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Abstract 



Financial time series typically exhibit strong fluctuations that cannot be 
described by a Gaussian distribution. In recent empirical studies of stock 
market indices it was examined whether the distribution P(r) of returns r(r) 
after some time r can be described by a (truncated) Levy-stable distribution 
L a (r) with some index < a < 2. While the Levy distribution cannot be 
expressed in a closed form, one can identify its parameters by testing the 
dependence of the central peak height on r as well as the power-law decay of 
the tails. In an earlier study [Mantegna and Stanley, Nature 376, 46 (1995)] 
it was found that the behavior of the central peak of P(r) for the Standard 
& Poor 500 index is consistent with the Levy distribution with a = 1.4. In 
a more recent study [Gopikrishnan et al., Phys. Rev. E 60, 5305 (1999)] it 
was found that the tails of P(r) exhibit a power-law decay with an exponent 
a = 3, thus deviating from the Levy distribution. In this paper we study the 
distribution of returns in a generic model that describes the dynamics of stock 
market indices. For the distributions P(r) generated by this model, we observe 
that the scaling of the central peak is consistent with a Levy distribution while 
the tails exhibit a power-law distribution with an exponent a > 2, namely 
beyond the range of Levy-stable distributions. Our results are in agreement 
with both empirical studies and reconcile the apparent disagreement between 
their results. 
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I. INTRODUCTION 



Financial time series are generated by complex dynamical processes that exhibit strong 
correlations between many degrees of freedom. The efforts to understand the dynamics of 
economic systems have involved empirical studies in which the temporal fluctuations of the 
prices of individual companies as well as of stock market indices such as the Standard & 
Poor 500 (S&P500) were examined These fluctuations can be characterized by the 

distribution of stock market returns as well as the volatility, that quantifies the magnitude 
of the market fluctuations. 

Consider a stock market index W(t). Its value is proportional to the average of the 
market values Wi, i = 1, . . . , N (given by the stock price of each firm times the number of 
its outstanding shares) of the N stocks that are included in this index. The fluctuations of 
W can be expressed in terms of the returns after a period of time r (say, in minutes), given 
by 

r(r) = lnW(t + r) - lnW(t). (1) 

For any value of r one can examine the distribution P(r) of the returns r(r). The number 
of independent data points available in the distribution is given by T/r, where T is the time 
period covered in the available data set. It was observed long ago that such distributions 
exhibit slowly decaying tails, unlike the Gaussian or exponential distributions. Moreover, 
the shape of the distribution was found to exhibit a self-similar form for different choices 
of t. It was proposed by Mandelbrot [[| that P(r) may be expressed by a Levy-stable 



distribution, L a (r), where < a < 2 | 10|,pJ| . Mathematically, the Levy distribution L a (r) 
is the limit n —>■ oo of the distribution of the sum of n independent stochastic variables 
taken from a power-law distribution of the form p(r) ~ r -1 ~ a when < a < 2 (that clearly 
exhibits an infinite variance). This is unlike the case of a distribution with a finite variance, 
that leads to a Gaussian distribution of the sum, according to the central limit theorem. The 
Levy distribution thus exhibits an infinite variance. However, in practical applications its 
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tail is truncated due to an upper cutoff in the power-law distribution that generated it |L2 . 
Although the Levy distribution cannot be expressed in a closed form |13| , it has two scaling 
properties that can be used in order to examine whether a distribution P(r) obtained from 
empirical data or numerical simulations is a (truncated) Levy distribution and to calculate 
its index < a < 2. The first property involves the dependence of the central peak height 
on the time r, that takes the form |13| 



L Q (r = 0)~r- 1 / a . (2) 

Thus, if the distribution of returns P(r) is a (truncated) Levy distribution, the value of a 
can be obtained from the slope of the graph of P(r = 0) vs. r on a log-log scale. The second 



property involves the power-law decay of the tails of the distribution that follows |~3| 



L a (r) ~ r- 1 -*. (3) 

Therefore, if the distribution P(r) is a (truncated) Levy distribution, the value of a can also 
be obtained from the slope of the tail of P(r) vs. r on a log-log scale. Obviously, a Levy 
distribution should satisfy the scaling relations for both the central peak and the tail, with 
the same exponent a. 

The distribution P(r) of the returns r(r) for the S&P 500 stock market index was recently 
studied for a range of r values, using the data for the six-year period of 1984-89 [|lj. The 
scaling of the central peak height vs. r was examined within the range of 1 < r < 1000 
minutes, yielding a straight line in the log-log scale over three orders of magnitude, with 
a slope that corresponds to a = 1.4. It was thus concluded that P(r) takes the form of a 
truncated Levy distribution L a (r) with the index a = 1.4. More recently the data set was 
extended to cover a 13-year period (1984-96) and was examined using the scaling analysis 
of the tail of the distribution P(r) of the returns r(r) for r in the range between 1 minute 
and 4 days. |2[. It was found that the tail of P(r) vs. r, on a log- log scale exhibits a straight 
line domain, indicating a power- law dependence given by Eq. (|3[). However, the slope was 
found to be consistent with a in the range 2.5 < a < 3.5, where the precise value depends 



on details such as the value of r and the fitting procedure. Clearly, these values of a are well 
outside the Levy-stable range of < a < 2. Therefore, not only that the distribution P(r) 
is not a Levy distribution with a = 1.4 - it is not a Levy distribution at all. Apparently, this 
result seems to be in disagreement with the conclusions of Ref. We thus observe that 
while the central peak maintains its Levy features the tails show a non-Levy behavior. In 
order to understand these puzzling results one needs to combine theoretical studies, suitable 
models and simulations of stock market dynamics, complementary to the empirical analysis. 

In this paper we study the distribution of the returns P(r) in a dynamical model that 
describes the time evolution of stock market indices [p^-|17|. The model consists of dynamic 
variables Wi, i = 1, . . . , N that represent the capitalization (total market values) of N firms. 
The dynamics represents the increase (or decrease) by a random factor \(t) [taken from 
a predefined distribution 11(A)] of the value Wi of the firm i between times t and t + 1. 
The dynamical rules also enforce a lower bound on the WiS, which is a certain fraction 
< c < 1 of the momentary average of the u^'s. This lower bound may represent the 
minimal requirements for a company stock to be publicly traded. It turns out that after 
some equilibration time the w^s exhibit a power-law distribution of the form p(w) ~ w^ 1 ^" 
[T7J. For any given value of N, the exponent a > is a monotonically increasing function 



of c. Since r(r) can be considered as a sum of r random variables taken from a power-law 
distribution p(w), one may expect it to converge to the Levy distribution L a (r) with the 
same exponent a. Since the power-law distribution is truncated from above, the tails of 
the resulting Levy distribution is also expected to be truncated [12]. Clearly, the dynamics 



is much more complicated. One reason for this is that the r random variables are not 
completely independent - they are taken from a finite set of N values of the w^s. Moreover, 
these values slowly change during the calculation of r(r), because at each time step one of 
the w^s is updated. 

To analyze the distribution of returns P(r) we first tune the parameter c (for the given 
value of N) to adjust the power-law distribution to the economically relevant case of a = 1.4 
jljl8l. We then examine the distribution of returns P(r) for a range of time intervals r and 



test the scaling behavior of the central peak as well as of the tails. It is found that the 
scaling of the central peak is consistent with a truncated Levy distribution with a = 1.4 for 
a broad range of 1 < r < 1000. For small values of r, up to about r = 50 (for N = 1000) the 
power-law decay of the tail of P(r) is also consistent with a truncated Levy distribution with 
the same value of a. However, for larger values of r the tail of P(r) exhibits a power-law 
decay consistent with a > 2, and thus deviates from the Levy distribution. These results 
are in agreement with the empirical analysis of the central peak presented in Ref. [|lj as well 
as with the more recent analysis of the tails presented in Ref. f2|. They thus reconcile the 
apparent disagreement between these two empirical studies. 

The paper is organized as follows. In Sec. II we present the model. Simulations and 
results are reported in Sec. Ill, followed by a summary in Sec. IV. 



II. THE MODEL 



The model [pl4|,|IT)|.|T7f describes the evolution in discrete time of N dynamic variables 



Wi(t), i = 1,...,N. At each time step t, an integer i is chosen randomly in the range 
1 < i < N, which is the index of the dynamic variable Wi to be updated at that time step. 
A random multiplicative factor X(t) is then drawn from a given distribution 11(A), which is 
independent of i and t and satisfies J x U(X)d\ = 1. This can be, for example, a uniform 
distribution in the range \ m in < X < X ma x, where \ m in and X ma x are predefined limits. The 
system is then updated according to the following stochastic time evolution equation 

Wi(t + 1) = X(t) Wi (t) 

Wj (t + l) = Wj (t) } j = l,...,N; j^i. (4) 

This is an asynchronous update mechanism. The average value of the system components 
at time t is given by 

1 N 



6 



The term on the right hand side of Eq. (|4]) describes the effect of auto-catalysis at the 
individual level. In addition to the update rule of Eq. (f|), the value of the updated variable 
Wi(t + 1) is constrained to be larger or equal to some lower bound which is proportional to 
the momentary average value of the wis according to 

Wi(t + 1) > c- w(t) (6) 

where < c < 1 is a constant factor. This constraint is imposed immediately after step (f|) 
by setting 

Wi(t + 1) — > max{w,(t + 1), c ■ w(t)}, (7) 

where w(t), evaluated just before the application of Eq. (f|), is used. This constraint describes 
the effect of auto-catalysis at the community level. Numerical simulations of the stochastic 
multiplicative process described by Eqs. (^) and (0), show that the w^s follow a power-law 
distribution of the form 

p(w) = Kw- 1 -" (8) 



for a wide range of lower bounds c, where K is a normalization factor [ 17 ] . It was found that 
the exponent a depends on the parameters c and N and is insensitive to the shape of the 
probability distribution 11(A). For simplicity, we use A uniformly distributed in the range 
0.9 < A < 1.1. 

III. SIMULATIONS AND RESULTS 

In the simulations below the number of dynamical variables is N = 1000 and the lower 
cutoff is chosen as c = 0.3, the value that provides the economically relevant distribution 
characterized by a = 1.4 [1,0. Under these conditions p(w) exhibits a power law distri- 



bution within three decades, between w min = 0.0003 and w max = 0.3. The data for this 
distribution was obtained from a large number of simulations collecting data at different 
times within each simulation after some equilibration time. To remove the possible effect of 



inflation, the values of the w^s fed into the distribution p(w) were normalized such that at 
any time t the sum J2i w i{t) = 1> namely w(t) = l/N [fT/j] . In the analysis of the returns, 
there is no need for such normalization adjustment, due to the fact that the returns quantify 
changes relative to the current value of w, namely they are normalized by definition. 

Consider the time evolution of the average w(t). At each time step, when Eq. (f|) is 
applied, neglecting the effect of the lower cutoff we obtain 

w(t + l)=w(t) + ±[\(t)-l} Wi (t). (9) 

This can be considered as a generalized random walk with step sizes distributed according 
to Eq. (H). Therefore, the returns after r time steps, given by 

r(r) = \nw(t + r) - \nw(t) (10) 



are expected [|12|] to follow a a truncated Levy distribution L a (r). Note that for small time 



intervals, the returns given by fliCf ) coincide with the relative change given by 

r(r) = " (t + r / ) 7" (t) . (11) 

However, for large r these two expressions provide significantly different results. 

In Fig. [l| we show the rescaled distribution r^P^/V 1 /") of the returns r(r) for r = 1, 
50, 200 and 1000. Near the central peak the four rescaled graphs collapse into a similar 
shape. The graphs for r = 1 and 50 maintain a similar rescaled form also in the tails while 
for larger values of r the tails go down more sharply. 

The value of a that characterizes the distribution can be obtained from the scaling of 
the central peak height as a function of r, according to Eq. (0). In Fig. |2] we show the height 
of the peak P(r = 0) as a function of r on a log-log scale. It is found that the slope of the 
fit is —0.71, which following the scaling relation of Eq. (0) means that the index of the Levy 
distribution is a = — 1/ (—0.71) = 1.4. 

To characterize the nature of the distribution P(r) we also examine the scaling behavior 
of the tails. For the Levy distribution the tail is expected to follow a power-law behavior 
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given by Eq. (|3|). In Fig. || we present the tail of the distribution P(r), on a log-log scale 
for t = 1. It is found that the slope is — (1 + a) = —2.4 which corresponds to a Levy 
distribution with a = 1.4. For larger values of r, the tails exhibit steeper slopes that exceed 
the domain of the Levy distribution, namely a becomes larger than 2. As an example, we 
present in Fig. ^ the distribution P(r) of r(r) for r = 10 4 on a log-log scale. We identify 
a range of about one order of magnitude in which the apparent slope is — (1 + a) = —3.5, 
namely corresponds to a = 2.5, which is outside the domain of the Levy distribution. It is 
thus observed that the tails of the distribution P(r) are much more sensitive to deviations 
from a Levy-stable process than the central peak. 

These results are in agreement with the empirical analysis of the central peak presented 
in Ref . |J as well as with the analysis of the tails presented in Ref . . They thus reconcile 
the apparent disagreement between these two empirical studies. To relate the parameters of 
the model more closely with the empirical studies we note that the typical time required for a 
single stock-market transaction is of the order of one minute. However, the transactions are 
done simultaneously in all the stocks included in the index that is analyzed. Therefore, the 
single transaction-time unit (say, one minute) roughly corresponds, in the model, to r = N 
time steps. The results of Fig. |41 for r = 10 4 are thus expected to correspond to a time 
interval of several minutes in the empirical analysis. Indeed, the value of a = 2.5 obtained 
in the numerical simulations is only slightly lower than the empirical results obtained for r 
in the range between 1 and 512 minutes. 

In the model we observe significant deviations from the Levy distribution as r increases 
towards the order of N. A possible explanation is that at this stage some of the W{S are 
already sampled more than once in a sequence of r time steps required to calculate one 
instance of r(r). This violates the requirement in the construction of a (truncated) Levy- 
stable distribution, that the r random variables should be independent. This starts to 
introduce significant correlations between the different variables that compose r(r). 

Another correlation effect is intrinsic to the calculation of the returns. Consider the 
return r(r), which is given by 
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r(r) — /Jn 



t=i 



(12) 



where the variable u^t) is independently picked at any time t. Note that the return depends 
on the normalized quantities w[ = Wi(t)/w(t). It is easy to see that the w'^s are not indepen- 
dent since at any time t they satisfy J2i w [ — N. This dependence is particularly apparent 
for the large w^s, since if one of them turns out to be extremely large the normalization 
condition prevents other w'^s from having values in its vicinity 



IV. SUMMARY 

Recent empirical studies of the fluctuations in stock market indices have provided con- 
flicting results. In these studies the distribution P(r) of stock market returns r(r) after time 
r were examined. The scaling of the central peak of P(r) was found to be consistent with 
a (truncated) Levy-stable distribution with index a = 1.4 |. However, the scaling of the 
tails, for a broad range of r values between 1 minute and a few days, was found to exhibit 
a power-law behavior with an exponent a = 3, which is well outside the range of the Levy 
distribution 0. 

In this paper we have examined the distribution P(r) for a model that describes the dy- 
namics of stock market indices. The model consists of dynamical variables u>j, i — 1, . . . , N, 
that describe the time-dependent market values of N firms, while their average is the corre- 
sponding stock market index. It was found that the scaling of the central peak is consistent 
with a Levy distribution and its index can be tuned to the economically relevant value of 
a = 1.4 by tuning a parameter. The tails of the distributions P(r) of the returns r(r), for 
a range of r values that corresponds to the empirically studied time intervals, were found 
to exhibit a domain of power-law behavior with a > 2, that falls outside the range of the 
Levy distribution. These results are fully consistent with the empirical results both for the 
central peak and for the tails and reconciles the apparent disagreement between them. 
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FIGURES 

FIG. 1. The rescaled distribution of the returns r 1 /«p( r (r)/r 1 / a ) for r = 1, 50, 200 and 1000. 
In the vicinity of the central peak we observe a collapse of all four graphs into a similar shape. 
The tails for the two smaller values of r follow the Levy-stable distribution with a = 1.4. The 
tails for the two larger values of r fall off more sharply and exhibit significant deviations from the 
Levy-stable shape. 

FIG. 2. The height of the central peak P(r(r) = 0) vs. r on a log-log scale. For a broad range 
of nearly three orders of magnitude in r values up to r = 1000, the slope of the straight line is 
— 1/a = —0.71, which corresponds to a Levy distribution with a = 1.4. 

FIG. 3. The distribution P(r) of r(r) on a log-log scale, for r = 1. The tail exhibits a range of 
power-law behavior according to Eq. (|3|) with a = 1.4, namely following a Levy distribution with 
the same value of a. 

FIG. 4. The distribution P(r) of r(r) on a log-log scale, for r = 10 4 . The tail exhibits a range 
of about one order of magnitude with an apparent power-law behavior. The slope in this range is 
consistent with Eq. (|3|) with a = 2.5. This value is not only different from the a = 1.4 observed for 
short times, but is outside the range for Levy-stable distributions. This curve strongly resembles 
the empirical distributions for the S&P500 presented in Ref. [2]. 
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